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ABSTRACT 

The phase diagram of the Schwinger model on the lattice with one and two degenerate 
flavours of Wilson fermions is investigated in the Hartree-Fock approximation. In case of a 
single flavour (not directly amenable to numerical simulation), the calculation indicates the 
existence of the parity violating phase at both weak and intermediate-to-strong couplings. 
In the broken phase, the Hartree-Fock vacuum sustains a nonzero electric fleld. With two 
flavours, parity is not broken at weak coupling. However, both parity and flavour become 
spontaneously broken at the Hartree-Fock level as the coupling becomes strong. 
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I. Introduction 



Understanding lattice fermions has been an outstanding issue since the beginnings 
of lattice field theory but a satisfactory insight is still missing. Aside from nontrivial 
numerical complications introduced by fermionic degrees of freedom, there are well-known 
conceptual obstacles, usually referred to as the "fermion doubling problem" . As revealed 
by the Nielsen-Ninomiya theorem [1], the doubling problem is intimately connected to 
chiral symmetry. In fact, this "No- Go" theorem has succeeded so far in preventing us from 
formulating chiral gauge theory on the lattice, although the investigations have intensified 
recently (see e.g. [2,3] and references therein) and possible clues might be at hand. 

Except from fundamental importance of chiral symmetry in electroweak theory, chiral 
symmetry has long been believed crucial in understanding the low energy behaviour of 
strong interactions, described by a vectorlike theory like QCD. Hinted by the small pion 
masses, the basic starting point is that the approximate chiral symmetry crucially shapes 
the way the low energy strong interacting world looks. Pions are regarded as Goldstone 
bosons coming from the spontaneous breakdown of this chiral symmetry, and the powerful 
predictions of the current algebra follow [4] . 

It is desirable to study these interesting issues within the nonperturbative framework 
of lattice QCD. Nielsen-Ninomiya theorem doesn't directly prevent us from "latticizing" 
a vectorlike theory, even in the chiral limit. However, the existence of the chiral anomaly 
complicates the situation considerably. The standard argument is that since the lattice is a 
physical regulator, any symmetry of the lattice action will remain the valid symmetry of the 
theory at every stage, including the continuum limit. Consequently, an explicitly chirally 
symmetric lattice model can not reproduce the anomaly structure of the vectorlike theory. 
Indeed, the most extensively used versions of lattice QCD, namely with Wilson fermions 
and Kogut-Susskind fermions, both explicitly violate chiral symmetry. Nevertheless, it 
can be shown in lattice perturbation theory that the amount of violation is just right to 
obtain the correct anomaly in the continuum limit [5]. This is of course quite comforting. 
However, to understand chiral symmetry on the lattice at finite lattice spacing (where all 
the numerical simulations are performed) means to relate the concepts of chiral symmetry 
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breaking to those of the lattice system, where chiral symmetry is not present. 

Restricting myself now to the case of Wilson fermions, it is generally believed that 
there is a line of phase transitions Kc{g) (in the hopping parameter - gauge coupling plane) 
running up from the QCD fixed point Kc{0)- On this line, pion-like state becomes massless 
and the continuum chiral QCD is believed to be approached by following this line towards 
Kc(0). Some time ago Aoki [6] set out to answer the questions raised in the previous 
paragraph in this context. Namely, since at ^ 7^ the masslessness of the pion can't be 
due to the spontaneous chiral symmetry breaking, what is it due to? In other words, what 
is the nature of the phase transition along K,c{g)'? 

Aoki's answer was that K,c{g) represents the line of phase transitions at which parity 
(one-flavour case) or parity and flavour (multi-flavour case) becomes spontaneously broken. 
For one flavour, the pion is identifled with the massless particle driving the parity violating 
phase transition {{iip'^Bip) 7^ 0). In case of two flavours, it is (i'i/'75'?'3'0) acquiring an 
expectation value, tt^ is identifled with a massless mode of this phase transition, while the 
charged pions are viewed as the Goldstone bosons coming from the breakdown of flavour. 
(■^75 IV') is always assumed to be zero thus giving no reason for 77 to be light, which could 
be viewed as a solution of the U(l) problem on a lattice [6]. 

At strong coupling, the supporting evidence for this scenario is quite convincing [6] . On 
the other hand in the weak coupling regime (relevant for the approach to the continuum), 
the situation is far from conclusive. Despite the fact that some numerical work has been [7] 
and continues to be done [8] , the very existence of the parity-flavour violating phase still 
needs to be examined, let alone the detailed picture of symmetry breaking. The existence 
of the parity-flavour violating phase has been recently established in the Nambu-Jona- 
Lasinio model using large N methods [9] and also in a numerical simulations with flnite 
number of colours (N=2) [10]. In that case, the parity-flavour breaking phase in a model 
with two flavours seems only to exist at strong and intermediate couplings. 

Recently, there has been a line of seemingly unrelated developments taking place 
concerning chiral symmetry following the ideas of Kaplan [11]. His approach amounts 
to the use of the surface modes of the vectorlike theory with Wilson fermions in 2d + 1 
dimensions as a basis for the construction of 2d-dimensional chiral gauge theory on the 
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lattice. (For review, see [12] and references therein.) Detailed Hamiltonian analysis of the 
surface modes of Wilson fermions has been carried out in [13], where it was also suggested 
that the notion of surface modes might be useful for understanding the conventional Wilson 
formulation as well. In particular, one can think of the appearance of the surface modes 
as an underlying mechanism generating the parity violating phase. 

The nature of the argument is as follows. Consider the Hamiltonian for one flavour 
of free Wilson fermions in one spatial dimension on a finite lattice with open boundaries 

Hw = K 5Z[^,+i(z7i - - V^,(i7i + r)V,+i] + M^Jj^j- (1) 
j j 

Here K, M, r are hopping parameter, mass and the Wilson parameter respectively. 

is a two component spinor living on site j. If one increases the hopping parameter (or 

equivalently decreases the mass) to the supercrticilal values, so that 



M 



2Kr 



< 1, (2) 



two levels start to behave differently from the rest of the spectrum and appear bound to 
the ends of the lattice [13]. As the size L of the system goes to infinity, the energy of 
these surface modes tends to zero. On a finite lattice the two modes mix and acquire the 
energy e ~ . Consider now the Dirac vacuum with all the negative energy levels filled. 
The last filled level will either be the surface mode on the right or the one on the left *, 
thus creating an asymmetric distribution of particles in the vacuum with an extra particle 
on one end. Consequently, after turning on the U{1) gauge field, this will generate an 
electric field running through the vacuum. In one spatial domension such a field can not 
be canceled by a pair production and we find that parity is not respected by the vacuum 
of this theory. 

These ideas were further developed in a recent inspirative review by Creutz [2]. He 
gives a comprehensive qualitative picture of the phase structure using the surface modes 
scenario in both single and multi-fiavour case. The argument is based on the frequently 



* On a finite lattice I assume that some very small left-right symmetry breaking term 
is added to the Hw, so that the left and right modes do not mix. 
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used analogy between nonabelian gauge theories in four dimensions and electrodynamics in 
two dimensions {QED2, massive Schwinger model). As is well known, QED2 exhibits some 
of the most intriguing features ascribed to QCD, namely confinement, chiral symmetry 
breaking and the existence of the 9 parameter. As such it represents a popular toy model 
for QCD with the advantage that it can be analysed semiclassically through bosonisation 
[14]. Taking lessons from the continuum and combining them with the existence of the 
doublers and surface modes on the lattice, one naturally arrives to the conclusion that the 
physics of the parity violating phase corresponds to the 6 = ti case in the continuum. This 
is based on the considerations of the previous paragraph and the fact that in QED2 the 
9 parameter has a direct physical meaning as the background electric field. The phase 
diagrams of Ref. [2] represent the expected positions of ^ = tt transitions on the lattice 
that should be applicable at weak coupling. 

The purpose of this paper is to investigate the issues of a parity violating phase 
in lattice QED2 with Wilson fermions in a direct lattice calculation. Compared to the 
wealth of exact and approximate information accumulated over the years on the continuum 
Schwinger model *, the knowledge we have on the lattice is rather modest. First of all, there 
are no exact solutions within any of the formulations where doublers are removed. The 
model was mostly investigated with Kogut-Susskind fermions, testing various numerical 
methods by comparing the lattice results to the exactly known continuum quantities in the 
massless case (see e.g. [16]). However, very little is known about the theory in the Wilson 
fermion formulation. The situation is particularly interesting for a single flavour, because 
in that case the direct numerical simulation is not possible. This is due to the fact that for 
certain gauge configurations the fermionic determinant is not positive and therefore one 
does not have a probabilistic weight for the purposes of Monte Carlo simulation. In a recent 
work [17], Gausterer and Lang studied the Lee Yang zeros of the partition function on small 
lattices analytically (at infinite coupling) and numerically (at intermediate couplings). The 
system was also studied in Ref. [18]. In case of two flavours, the direct numerical analysis 
is possible, but to my knowledge, the systematic study of phase structure has not been 
carried out. 

* The recent activity concentrated mostly on the multi-flavour case. See e.g. [15]. 
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In what follows, I will study the phase diagram of Hamiltonian QED2 on the lat- 
tice with one and two degenerate flavours of Wilson fermions in the Hartree-Fock (H-F) 
approximation. I will work in the axial gauge, where the gauge degrees of freedom are 
easily eliminated in favour of fermionic fields, thus providing a convenient setup for the 
use of the H-F approximation to the vacuum of the theory. The H-F ground state, or the 
independent fermion ground state, is a state one might call (by definition) the mean field 
ground state for the theory of interacting, particle number conserving (charge conserving) 
fermions. It has a nice variational interpretation and can be regarded as a first term in 
a series of systematic variational improvements. While belonging to the standard set of 
techniques used in many-body theory, H-F methods are rarely invoked in the lattice gauge 
context. In the case under consideration, however, it can give us valuable hints about the 
phase structure of the theory. 

Investigating the phase diagram, I adopt the point of view, taken in Ref. [2] , that it is 
quite natural to introduce the axial mass term (the "M5"-term). More specifically, I will 
consider the generalized version of Hamiltonian (1), given by 

Hw, =Hw + M^i^j^^i/jj = Hw + H5. (3) 

Indeed, one obtains such a term from the conventional mass term by a chiral rotation. 
The existence of a ^ parameter in the continuum theory can be thought of as being due 
to the fact that because of regularization, the theory with such a rotated mass term is 
actually not equivalent to the original one and is thus anomalous. In this way, introducing 
M5 essentially means trading 6 (being an independent parameter of the theory) in favour 
of this new mass term. Consequently, I will consider the lattice theory in the space of 
three bare parameters, M, M5 and g (gauge coupling), in contrast to three conventional 
continuum parameters m, $ and e. 

In Sec. II, starting from the continuum theory in axial gauge, I will formulate the 
lattice model with one flavour and discuss its discrete symmetries. The H-F approximation 
is then described in Sec. III. I do this in some detail, stressing the variational character of 
the method and the fact that all the discrete symmetries (including parity) are preserved 
by the approximation. The numerical solutions of the H-F equations on flnite lattices and 
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their implications on the H-F phase diagram are discussed in Sec. IV. I try to make all the 
observed qualitative features of the vacuum plausible by tracing them to the elementary 
picture of interacting particles in filled H-F levels. The model with two flavours is then 
formulated and analysed in Sees. V and VI. Summary, together with some generalizations 
and speculations, is given in Sec. VII. Finally, the numerical procedure used to solve the 
H-F equations is described in Appendix. 



II. The Model with One Flavour 



The massive Schwinger model with one flavour of fermions is defined by the Lagrangian 
C = —F^.F^- +^r{id^^ - eA^) - m] , (4) 

where F^i, = dfj,Ai, — di,A^. ij; is the two-component spinor field, is the gauge field and 
m, e the mass and dimensionful coupling constant respectively. I will consider the theory 
in axial gauge {Ai = 0). 

The equation of motion for Aq is the equation of constraint (Gauss' law) 

E'{x) = -A'^{x) = eij\x)ij{x) = ep{x) , (5) 

where E{x) = —Aq{x) is the electric field and p the charge density. It can be solved by 

Ao{x) — —- J dx' I X — X I ep(x') + Bx + C . (6) 

The integration constant C is physically irrelevant and will be put to zero. The constant 
B represents a uniform background field which can be also put to zero for my current 
purposes. Then the Hamiltonian of the theory for the zero total charge takes the form 

,2 



H = J dx V'l^Ti^i + 1Tl)^jJ — —J dx'p{x) \ x — x' \ p{x') 



(7) 



Upon quantization, if) and H become operators in the corresponding Hilbert space, 
with field operators subject to the canonical anticommutation relations. In the standard 
treatment, the local charge density operator is replaced by its normal ordered version 



: tp'^tp:, with the normal ordering usually performed with respect to the filled Dirac sea (no 
particles, no antiparticles) . This effectively amounts to the compensation of the infinite 
charge, generated by the sea. Since I will not introduce the antiparticlc operators it is 
more convenient for my later purposes to define the charge density in such a way that the 
compensation is explicit, namely 

p{x) = ip^{x)il;{x)-l. (8) 

To formulate this theory on a lattice is now straightforward. It will be defined by the 
Hamiltonian 

H = - Pn \ n - m \ pm = Hw^ + Hi , (9) 

n,m 

where is the free part with Hr, defined in (3) and 

Pn^^i^n-l- (10) 

The indices n, m label the lattice sites runing from 1 to L and the lattice spacing has been 
set to unity. Fermionic variables are subject to open boundary conditions and satisfy the 
canonical anticommutation relations 

{V'-,Vj:} = (11) 

with cr, T being the spinor indices. The electric field in this formulation is a derived quantity, 
defined through the the lattice analog of Gauss' law (5) by 

^.- = f[E«- E«]- (12) 

1=1 1=3 + 1 

Here Ej is the operator of electric field on link (j, j + 1). 

It is worth emphasizing at this point that similarly to the continuum case, the lattice 
interacting theory posesses exact discrete symmetries. Choosing the representation of 7 
matrices as 

V=(_°, J) .. = .V, (13) 
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the operations of C,P,T in continuum (left column) and on the lattice (right column) are 
defined by the following transformation properties of the field operators 

Ci(;{x)C-^ = -f{f{x) C'^jC-^ = 7iV^J (14) 

Til;{x)T-^ = joi^{x) Ti^jT-^ = 7oV'i • 

Here C, P are unitary and T antiunitary operators. The above transformations indeed leave 
the corresponding continuum and lattice Hamiltonians without (M5 = 0) unchanged. 

Note that H5 is odd under P, C and even under T. Therefore CP and T are the 
exact symmetries of the theory at any M5. This is different from the situation in four 
dimensions, where is invariant under C and changes sign under P, T, which in turn 
implies that both PC and T are explicitly broken by nonzero M5. 

III. Hartree-Fock Approximation 

In solid state physics, the H-F approximation is frequently referred to as the "inde- 
pendent electron approximation" and this probably captures its essence best. Indeed, the 
main idea is to approximate a given state of the fermionic many-body system by a Slater 
determinant of some set of one-particle states. This is usually applied to approximate the 
unknown vacuum of the theory, which is also my main interest here. In that case the 
method boils down to finding a set of one-particle states with the Slater determinant of 
minimal energy. 

Let me therefore start to investigate the vacuum of the one-flavour lattice Schwinger 
model by considering an arbitrary (but fixed) complete orthonormal set of one-particle 
fermionic states on a lattice of L sites 

S' = {0" I a = 1,2,...,2L}. (15) 

Every state (j)'^ is a collection of two-component spinors, residing on site n 

r^iK I n = l,2,...,L} (^S)- (16) 
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One can build the fermionic many-body Fock space out of these states in a standard way 
and define the complete set of fermionic annihilation operators by 

V'n = 5]aaC • (17) 

a 

With definition (10) of the local charge operator, restriction to the charge zero sector 
translates into the requirement of half-filling. In other words, only states with L particles 
are allowed. An arbitrary state | if^) from this subspace can be written in the form 

I V') = c° I 0/,07) + Y^cl I 1/,17), + Y^cl I 2/,27), + ... + I L/,L7). (18) 

r r 

Here | 0/, 0/ ) , the "sea" , is the state with first L levels filled 

|0/,07) = al4...a[|0). (19) 

I n/, nf ) generically represents states with n fermions removed from the sea and put into 
n empty levels (n fermion - n antifermion states). Index r enumerates these basis states 
at fixed n. Obviously, the dimension of the charge zero sector is (^^) . 

Apart from charge conservation, no other symmetries are assumed to be respected by 
the vacuum. One usually restricts the space of states further by going to subspace of zero 
momentum. However, momentum is not well defined on a lattice with open boundaries. 
Specification of (^^) complex coefiicients in the above decomposition therefore consti- 
tutes an exact representation of the vacuum. Requirement of minimal energy defines a 
variational problem for determination of these coefficients. However, on a reasonably sized 
lattices the number of variables in the problem becomes too huge to be manageable. 

The variational philosophy behind the H-F approach is to regard the set S' as a col- 
lection of variational parameters instead. In particular, the method aims at adjusting the 
one-particle basis in such a way that while retaining only the first term in decomposition 
(18), the lowest energy state is achieved. Note that in this way, the variational problem 
involving 2(^^) real variables is replaced by one involving 3L^ real variables *. Moreover, 

* Counting here includes the fact that only the filled levels represent the true variational 
variables since only they contribute to the total energy. One also has to take into account 
the orthonormality constraints. 
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fixing S by the Hartree-Fock prescription transforms the decomposition (18) into a well 
defined variational improvement scheme. Indeed, employing the H-F basis one expects the 
states I n/, nf ) to play increasingly less important roles in the true vacuum with increasing 
n. This is expected to be true regardless of the value of the gauge coupling. 

Having defined the H-F approximation it is now a straightforward matter to transform 
the problem into the familiar manageable form. Using the field decomposition (17) one 
can rewrite the free part of the Hamiltonian (9) as 

Hw, = Y.'i'nKnm^m = J] /C"^a;3 , (20) 
nm al3 

where 

nm (21) 
Knm = K [5n,m+llo{ill - r) - 5n+l,m7o(^7l + ^)] + ^nm [M70 + iM^-^] . 

The interaction term becomes a little more complicated and has three parts 

Hi = ^i+ y E^a-^"''^^ + C 5^ aia^^M-^^'a,as , (22) 

a/3 a/375 

where 

■A^'^^ = E <^n^ I ^ - ™ I = E ^nUZ \n-m\ cf>^Jcf>i^ (23) 

nm nm 

2 

and £1 is an unimportant constant £1 = \ — 'm- \- Note that loosely speaking, 

the constant term corresponds to the self-interaction of the Dirac sea-compensating charge, 
the quadratic term arises due to the interaction of this charge with the system and the 
quartic term represents the interactions of the system itself. 

The mean energy in the Slater determinant (19) is a function of the set S and is given 

by 

L 2 2 ^ 

£{S) = £1 + Y, [''C"" + ^AA""] + ^ E [-^"^^^ - . (24) 

a=l cep = l 

The Hartree-Fock set of states S^^ is now determined by minimizing £{S) with variables 
0^ subject to the orthonormality constraints. Standard manipulations then reveal that 
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this variational problem can be solved by subjecting the one-particle wavefunctions to the 
H-F equations of the form 

Y^HnmCt^'^a = ^"C , (25) 



m 



with 



and 
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^ nm ' ^ nm 



(26) 



VZ. = Snm E I ^ - i I (1 - E ^T^^) ^nm = Y.<i>i\n-m \ cf>^J . (27) 
j f3=l /3=1 

The "twidle" in Hnm serves to denote the fact that these are not the one-particle matrix el- 
ements of Hamiltonian (9). As expected, the H-F equations take the form of a one-particle 
Hamiltonian eigenstate problem with the complication that the Hamiltonian matrix de- 
pends on the eigenstates themselves. Thus the equations have to be solved self-consistently. 
As usually, the self-consistent potential has direct and exchange parts. 

Due to the self-consistent feature of the above H-F equations, finding an exact solution 
is a nontrivial task and I haven't succeeded in doing that. On the other hand, there is a 
simple way to attempt to solve these equations on finite lattices numerically, namely by 
iteration. Straightforward application of the iterative procedure however doesn't converge 
to the self-consistent solution. The nature of the problems is similar to those described in 
Ref. [19] in the context of continuum QCD2 in th large N limit. I discuss these technical 
issues in Appendix. Using a modified approach, numerical solutions can be iteratively 
found in wide range of coupling constants. 

Let me close this section with a few remarks concerning symmetry within the H-F 
approximation. It is usually helpful and desirable that the approximation scheme retains 
as much symmetry of the approximated system as possible. Especially if the main purpose 
of the investigation is spontaneous symmetry breaking. In particular, as discussed in the 
previous section, if M5 = 0, the Hamiltonian (9) is invariant under parity. Is the parity 
invariance present in the approximation? In what sense? 

The underlying dynamics driving the H-F approximation is entirely embodied in the 
H-F equations. Therefore, the symmetries of these equations should also determine the 
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symmetries of the approximation. It follows from transformation properties (14) that 
under the operation of parity the one-particle wave function Xn transforms into joXL+i-n- 
Since the H-F Hamiltonian implicitly depends on its eigenstates, the parity operation in 
this case has to involve the whole set S. In fact, it can be checked quite easily that the 
operation 

•5^^ = m PS'''' = {70Cn} , 6" ^ 6" (28) 

is a symmetry of the H-F equations. In other words, if is a self-consistent set solving 
(25), then PS''^ is also a self-consistent set with corresponding one-particle energies equal. 

Let me also mention that not just parity but all the discrete symmetries discussed in 
the previous section are preserved by the H-F approximation in the above sense. Of course, 
performing the symmetry operation on the H-F vacuum can lead to a Slater determinant 
involving different self-consistent set, thus opening the possibility of spontaneous symmetry 
breaking. 

IV. Numerical Analysis (One Flavour) 

In this section, I will discuss the results of the H-F analysis for the model with one 
flavour. I will concentrate on the phase diagram in the M — g^ plane. To observe the parity 
violating effects on a finite lattice, I fix M5 to a very small value (M5 = 10""^) throughout 
this section. The values of hopping parameter {K — 1) and Wilson parameter (r = 0.5) 
are set in such a way that the critical value of M at zero coupling is Mc(0) = 1. All 
quantities are given in the lattice units. Moreover, the electric field is allways measured in 
units of g. In particular, this expectation is calculated using formula (12) with operators 
of charge density replaced by their expectation values in the H-F vacuum, 

(Pn)^^ = E<^nVn ' 1, (29) 

and the factor of g removed. 
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To assess the accuracy of the self-consistent numerical solution and to provide the 
criterion for terminating the iterative procedure, I have computed 

5^=) = ^ I (/.^C^+i) - 0«W I (30) 

a,n 

at each iterative step k. Here ^n^'^^ £ S^^\ the one-particle set after k iterations (see 
Appendix). Obviously, 5 = only for the self-consistent set. In all cases discussed in 
this section S^'^^ has been accepted as a numerical solution only if 6^'^^ < 10~^. In most 
cases however, this number has been much smaller (up to four orders of magnitude at 
weak couplings). With the above bound, the physical characteristics of the H-F Slater 
determinant (such as energy), became essentially insensitive to further decrease of 5. To 
achieve this accuracy on the lattices I have studied (L = 32, 40, 48) took typically a few tens 
of iterations at weak couplings {g'^ < 0.5) and a couple of hundreds at intermediate and 
strong couplings {g"^ > 0.5). Working in the vicinity of the phase transition typically added 
roughly one order of magnitude to the number of iterations. In the region of couplings 
studied here (^^ < 3), the self-consistent solution has always been straightforwardly found 
with free wave functions at given M used as a starting point for the iteration. 

The representative example of the most relevant finding in this study is displayed in 
Fig.l. The vacuum expectation value of the electric field in the middle of the 32-site lattice 
is plotted as a function of the fermion mass at weak coupling {g"^ = 0.1). Note that for 
large values of M, the electric field tends to zero as one would expect in the parity-invariant 
theory. However, at small fermion masses the field acquires an expectation value and the 
two regions are separated by a rapid transition. This suggests the existence of a parity 
violating phase transition and confirms the qualitative picture presented in Ref. [13] at the 
Hartree-Fock level. 

The spatial dependence of the electric field across the lattice is plotted in Fig. 2a for 
typical cases in the broken and symmetric phases. Note that in the broken phase, the 
field nicely settles to a uniform bulk value essentially across the whole lattice. In the 
symmetric example, the field is almost zero everywhere. It is quite interesting to see the 
spatial distribution of Hartree-Fock levels in these two situations. This is shown in Fig. 2b 
where I plot the energy of these levels against the mean position of particles in them. In 
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the symmetric case, all the particles reside on average in the middle of the lattice and the 
left-right symmetry is preserved up to small explicit violations caused by the presence of 
the small Ms-term. The filled levels produce a uniform charge distribution neutralized by 
the compensating charge. This is to be compared to the situation in the broken phase, 
where the left-right symmetry is completely lost. Indeed, it is energetically favourable for 
the levels to spread out asymmetrically. Filling the sea generates the surface charge and 
an electric field. 

Similar behaviour is observed also at higher values of the gauge coupling. The resulting 
positions of phase transitions observed on the lattice with 48 sites are plotted in Fig.3. 
The transition points here are determined simply as the locations of the rapid rise of the 
vacuum expectation value of the electric field at fixed coupling. In particular, the phase 
transition is assumed to happen at the fermion mass Mdg"^), where this expectation value 
rises above 10~^, i.e. above the value one order of magnitude larger than the size of the 
parity violating Ms-term. By comparing to the results on smaller lattices (L = 32, 40) I 
expect the critical masses at nonzero couplings to be increased by a few parts per hundred 
in the infinite volume limit. 

For the model in the standard Euclidean formulation, Gausterer and Lang [17] con- 
cluded the existence of a phase transition at infinite coupling. After appropriate rescaling 
the parameters of their model, the quoted position of this transition is Mc(oo) ~ 0.32. 
Although I don't know of any apriori reason why the phase transition should occur at 
the same place in both formulations, it is interesting to observe that their result is an 
acceptable asymptotic value at strong coupling here too. 

Similarly to the electric field, the simplest local fermionic parity-odd operator, namely 
axial charge density ii/j^^i/j, also acquires an expectation value at the parity violating phase 
transition. This is illustrated in Fig.4a where I plot both the electric field and the axial 
charge density as a function of fermion mass at g'^ = 1.0. Both operators appear to acquire 
an expectation value simultaneously as expected. Typical spatial dependence of the axial 
charge density in the broken and symmetric phases is plotted in Fig. 4b, showing the bulk 
nature of the order parameter. 

The relative size of the electric field and the axial charge density in the broken phase 

15 



varies with gauge coupling. This is demonstrated in Fig. 5 where I plot these expectation 
values at fixed fermion mass. Note that while the electric field starts up finite at weak 
coupling and decreases monotonically at intermediate and strong couplings, axial charge 
density behaves in a complementary way. It approaches zero with vanishing coupling and 
rises as the coupling increases. 

The above behaviour of electric field is simply a manifestation of charge shielding, an 
effect well known to be present in the continuum theory as well. Indeed, consider first 
the free theory. In that case it is just the filled surface mode that is responsible for the 
parity-breaking effects. That's why the electric field approaches value 0.5 {9 = tt) when 
coupling tends to zero. However, once the gauge coupling is turned on, the rest of the 
levels spread out (see Fig. 2b) and the accumulation of surface charge is a result of the 
collective action of all self-consistently interacting particles in filled states. The net effect 
of this phenomenon is a screening of the surface charge. As the value of the gauge coupling 
increases, while remaining in the broken phase, one expects the levels to spread out and 
screen even more since the system wants to reduce the positive attraction energy of the 
surface charges. For example, at g'^ = 1.5, the spatial distribution of H-F levels is shown 
in Fig. 6. At strictly strong coupling, when the interaction term absolutely dominates, 
every particle in the sea will live bound to just one site of the lattice, thus eliminating the 
surface charge completely. Therefore, the field is expected to vanish in this limit in the 
H-F approximation. 

In the light of the above considerations, behaviour of the axial charge density in 
Fig.5 becomes also quite natural. Indeed, since at zero coupling the parity violation is 
all concentrated on the ends, it will not be refiected in the expectation value of the local 
operator inside the system. Consequently, one expects the bulk axial density to vanish. 
At strong coupling however, parity violation is equally contributed by all the filled levels 
and the axial density acquires an expectation value. 

An interesting feature already present in the examples of Fig. 2, but quite striking in 
Fig. 6, is that because of the interaction energy, it is not necessarily the lowest one-particle 
states that are filled to form the H-F vacuum. Indeed, in Fig. 6 almost half of the filled levels 
(denoted by diamonds) are those with positive one-particle energies. This is discussed in 
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more detail in the Appendix. 

Finally, let me close this section by discussing the order of the parity violating phase 
transition in the H-F approximation. It is well known that the mean field-like approaches 
are frequently misleading about the order of the phase transition and the critical exponents. 
Therefore, while I assume that it is plausible for H-F approximation to recognize the 
transition, the information it gives about the order should be taken with some care. 

I have calculated the connected correlation functions for electric field and axial charge 
density in the H-F vacua. For the case of electric field, these functions have a very nice 
exponential decay in all cases I have studied and the corresponding correlation lengths 
could be reliably determined. Typical behaviour of the inverse correlation length (mass 
gap) across the phase transition is shown in Fig.6a. In the symmetric phase, the correlators 
of the axial charge behave in the same way. However, in the broken phase they show some 
differences as can be observed from Fig. 6a. For weak couplings, the determination of the 
correlation length from these axial charge correlators in the broken phase was less accurate 
than from the electric ones. Using the electric field correlation functions. Fig. 6b shows the 
mass gap along the line of phase transitions on a lattice of 48 sites. These results are 
reasonably finite-size stable (more so at stronger couplings) and I exclude the possibility 
of mass gap reducing to zero in the infinite volume limit. I conclude that in the H-F 
approximation the phase transition is of first order at finite approaching a second order 
endpoint at zero coupling. 

V. Two Flavours 

Turning now to the case of two degenerate fermion fiavours, I will consider the lattice 
Hamiltonian 




(31) 



/ nm 



fh n,m 



Here f,h are flavour indices assuming two values, K^m is defined in (21) and 



(32) 
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The fermionic operators V'n subject to the canonical anticommutation relations. 

Note that with flavours being degenerate, the parameter space of this two-flavour the- 
ory is the same as for a single flavour, namely M, M5, g. Also, similarly to the one-flavour 
case, the above lattice model retains all the discrete symmetries of the corresponding 
continuum theory. In addition, the two-flavour model is invariant under unitary transfor- 
mations in flavour space. All of these symmetries will be preserved by the Hartree-Fock 
approximation in the sense discussed in Sec. III. 

Inclusion of the flavour index does not require any conceptual changes in the ap- 
plication of the Hartree-Fock procedure. On the technical side, it is easiest to skip the 
explicit use of flavour notation and assemble the two fermionic operators on site n into a 
4— component column ■^^ = (^5) . The Hamiltonian then takes the form 



Knm=(^r J" ) (34) 



H = ^ IplKnm'ipm - — ^ pn \ n - m \ Pm , (33) 

nm n,m 

with 

— ( ^nm 

K 

and 

p„ = iplipn - 2 . (35) 

Note that to compensate for the charge of the Dirac sea now requires two units of charge 
per site. 

Using the above notation and the decomposition 

4L 
a=l 

the derivation of the H-F equations is a line by line repetition of the procedure for the 
one flavour case up to the factors of two coming from the doubling of the compensating 
charge. Indeed, the H-F equations take the form 

J2 Hnm^m = ^"C , (37) 
m 

with 



~ _ 9^ 

Hnm — ^nm ~l~ 
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' nm ' ' nm 



(38) 



and 

2L 2L 

K^'^ = 5^m$^|n-J I (2-5^(/>fV^^) V;^^ = ^</>^ |n-m|</>^t. (39) 

J ^=1 /3=1 

Note that these equations are formally almost identical to those for one flavour. The 
crucial diff'erence however is that (/)^ is now a 4— component object and H a 4L x 4L 
complex matrix. 

VI. Numerical Analysis (Two Flavours) 

I will start with the discussion of the phase structure in the M — M5 plane at fixed 
gauge coupling. Before turning to the results of the H-F analysis, let me first briefly explain 
what one would expect to be happening here at weak coupling based on the surface mode 
picture. To do that, I will need to borrow the of ideas of Ref. [2], and to make this paper 
reasonably self-contained, to review briefly the part that is relevant here. Thinking first 
in the continuum context, consider the standard mass term m'4>'4> and its transformation 
under the chiral rotation ij) — > e'^^^^'^^ip. We have 

mipip — )• mcos{9)'ip'ip + msm{9)i'ip^5'ip . (40) 

Therefore, the chiral rotation by angle 9 corresponds to the rotation of the vector (m, 0) 
in the m — m^ plane around the origin ("chiral point") by the same angle 9. While naively 
expecting that the physics should be the same after the above change of variables, this 
is actually not the case because of the chiral anomaly. What we are actually getting is a 
physics with different "gauge" ^-parameters, i.e. with different background electric field, 
realized for example through the existence of the surface charges. With this identification, 
the above transformation prescription gives the approximate relation (the renormalization 
effects, for example, will shift the chiral point to negative m) between the theory considered 
in the parameter spaces (e, m, 6) and (e, m, m^). 

On a lattice with Wilson fermions, the situation is a little more complicated, because 
except from (Mg, 0), there is another chiral point in the M — M5 plane, namely (— Mg, 0), 
where the doubler goes massless. The conjecture then is that here the total ^-parameter 

19 



gets two contributions, each being the angle with respect to the two chiral points, with 
doubler's contribution taken with the reversed sign. These angles are sketched in Fig. 8a. 
Nf degenerate flavours will contribute equally to the total value of 9 and consequently, 
one expects the following approximate relation to hold at weak coupling 

e = Nf{ep-ed). (4i) 

Here 9p, 9 a are the contributions of the particle and the doubler respectively. In other 
words, on a lattice with open boundaries, the system is expected to generate surface 
charges in such a way, that the resulting electric field will approximately correspond to 9 
given by the above relation. 

If this qualitative picture is correct, there should be phase transitions occuring in 
the M — M5 plane at the positions where 9 reaches n. Indeed, for ^ > tt, it will be 
energetically favourable to create a fermion-antifermion pair thus reducing the magnitude 
of the electric field and switching its sign. Consequently, 9 should jump from tt to — tt 
across these phase transitions. The condition 9 = tt defines a line in M — M5 plane, but 
its qualitative behaviour strongly depends on the number of flavours. With single flavour, 
the only solution is a straight line, connecting and —Mc- This can be understood 
already from the point of view of the "naive" surface mode picture as I presented in the 
Introduction. Indeed, switching the sign of a small parity violating "M5"-term causes the 
two surface modes to exchange the ends of the lattice, thus switching the sign of the surface 
charges and the electric fleld. For two flavours, the condition 9 = tt deflnes a circle with 
centre at the origin and radius Mc- Including a small flavour breaking in both M and M5 
to visualise the chiral endpoints, the phase diagram is expected to look qualitatively as 
sketched in Fig. 8b [2]. Note that contrary to the single flavour case, the prediction of this 
phase structure in the M — M5 plane is quite nontrivial. 

Guided by this simple picture at weak coupling, I set out to look for these phase 
transitions in the Hartrcc-Fock approximation. For numerical work, I again flxed the 
values of the hopping parameter (K = 1) and the Wilson parameter (r = 0.5), so that the 
critical value of M at zero coupling is Mc(0) = 1. Also, I have always included a very small 
explicit flavour breaking in both M and M5. In particular, the masses of flavours were 
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of the form M ± 0.001 and M5 ± 0.001. The self-consistent solutfons were again obtained 
by the modified iteration procedure as described in Appendix. For all results presented 
here the accuracy of the solution, given by Eq. (30), was better than 10""^ and typically 
about 10~^. To achieve this accuracy in the immediate vicinity of the phase transition on a 
lattice of 40 sites took less then 10^ iterations at strongest coupling studied here {g'^ = 6). 

The graphs in Figs.9a,b illustrate how the above qualitative ideas are reflected in the 
Hartree-Fock approximation. In Fig. 9a I plot the vacuum expectation value of the electric 
field along the Ms-axis (M — 0) at = 0.1 on a, lattice with 40 sites. Note that the H-F 
vacuum nicely exhibits the expected abrupt change in the electric field and the reversal of 
its sign. Probing the field along the M-axis (with small M5 present) gives the dependence 
ploted in Fig.Qb. While behaving in qualitatively the same way, the magnitude of the field 
is becoming small as one approaches the M — g"^ plane. This is what one would expect 
if these 9 = tt transitions, which are naturally first order, end in a second order chiral 
endpoint in this plane. The anlysis of the electric field correlators suggests however, that 
in the H-F approximation, these transitions at finite coupling are first order even close to 
the M — g'^ plane. It should be stressed again however, that this might well be an artifact 
of the approximation. 

Defining the transition point as the position in the M — M5 plane where the field 
switches its sign, I plot the phase diagrams for ^r^ = 0.1 and ^r^ = 1.2 on a lattice with 32 
sites in Fig. 10. In fact, only the points in the upper right quadrant were really calculated. 
The rest of them were obtained using symmetry with respect to the mass reflections. Note 
that while not exactly of circular shape, the transition lines reflect the expected qualitative 
features deduced from the surface mode picture. Also, at (7^ = 0.1, non- negligible finite- 
size effects are present here. While the transition point along the M-axis is essentially 
stable against the increase of the lattice size, the transitions along the Ms-axis occur at 
M5 = 0.36, 0.42, 0.48 on the lattices with 32, 40 and 48 sites respectively. Thus it is quite 
possible that the ellipse-like shape of the transition line will become more circle-like in the 
infinite volume limit. 

Let me now turn to the question of Aoki's phase in the M — g"^ plane. First note that 
the surface mode picture doesn't suggest that parity-fiavour broken phase should exist here 
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at weak coupling. Indeed, as one turns off the M5 and moves along the M-axis, the angle 
9 defined by (41) is always zero. Both flavours will generate their surface mode as M is 
lowered below Mc, thus changing 6 by 27r and physically changing nothing. A good way 
to picture this is by looking at Fig. 8b: However small the flavour breaking is, if one moves 
close enough to the M-axis, it is always possible to pass below the chiral point. 

Discussion of these issues in two dimensions might appear a little academic since 
because of the Mermin- Wagner theorem [20], one would not expect flavour to be spon- 
taneously broken here at any coupling. Nevertheless, if the surface mode picture is of 
relevance to QCD in four dimensions, then this qualitative prediction wouldn't change. 
Moreover, I consider the following to be good reasons to investigate these issues in the 
context of the Schwinger model itself: 1) The lattice model, deflned by (33), is nonlocal, 
and as such does not exactly satisfy the usual assumptions of the theorem. 2) Parity can 
still be broken. 3) Flavour could be erroneously broken within the H-F approximation. 
Then, although not useful as an information about the lattice Schwinger model, it can 
serve as a toy picture of what might be happening in QCD where there is quite convincing 
evidence that parity-flavour is broken at strong enough coupling. 

With that in mind, I have calculated the expectation values of '075 ''"3 and i/'75l'i/' in 
the H-F vacua. Here is the third Pauli matrix and 1 a unit matrix in flavour space. 
Note that the form of flavour breaking in M5 used here chooses the rs-direction if flavour 
is broken. Note also that if M5 = and (V'Ts'^sV') 7^ 0, it is both parity and flavour that 
are spontaneously broken. On the other hand, if this expectation is nonzero at nonzero 
M5, the parity is broken explicitly while the flavour spontaneously. Furthermore, if at 
M5 = we had ('075l'0) ^ and {4'75'!'3'4^) — 0, it would indicate that only parity has 
been spontaneously broken. 

The results of the H-F analysis in M — plane (with M5 = 10~^) on the lattices with 
up to 40 sites are as follows. I have found no evidence of V'TslV' acquiring an expectation 
value in the region of couplings < 6. Consequently, there is no indication of parity being 
broken alone. However, there are regions where {ip'j^Tstlj) is nonzero in the Hartree-Fock 
approximation on the flnite lattice. For example, in Fig. 11a I plot this expectation as a 
function of fermion mass at = 1.2. The broken region appears as a narrow peak adjacent 
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from the left to the = tt" * transition point on this lattice of 40 sites. I have observed 
similar peaks at g'^ = 0.1 and g'^ = 4.0 with heights roughly 0.08 and 0.70 respectively. 
For the two weaker of the above couplings, I have also performed a finite-size analysis of 
the width of this broken region. This width decreases linearly with exhibiting a small 
negative intercept in both cases. On the other hand, the heights of the peaks stay constant 
as the lattice size increases. I therefore conclude that these narrow regions will not survive 
in the H-F phase diagram in the infinite volume limit. The only remnant of them will 
probably be the singular behaviour of the parity-fiavour order parameter at the "6* = tt" 
transition point. 

The situation qualitatively changes at even stronger couplings. In particular, the 
parity- flavour broken phase indeed opens up at the subcritical fermion masses. This is 
demonstrated in Fig. lib, where I show the behaviour of the order parameter at g'^ = 6.0 
on a lattice with 32 sites. There has been a negligible change here as the lattice size 
increased to L = 40. I therefore expect the finite-size effects to be small. It is also worth 
mentioning that (V'Ts'T'sV') 7^ in the whole inside region of the "9 = tt" line at this strong 
coupling and not only in the M — g'^ plane. This is in contrast to the case of narrow broken 
regions at weaker couplings which can only be observed close to M — plane on a finite 
lattice. 

This concludes the review of the most important aspects of the numerical information 
obtained in this study. In the last section, I will turn to generalizations and speculations. 

VII. Summary, Generalizations and Speculations 

The Schwinger model on a lattice with Wilson fermions has been studied in the 
Hartree-Fock approximation. The main focus was given to the global structure of phase di- 
agrams with one and two degenerate flavours of fermions. In future communication, I plan 

* Note that I loosely refer to a transition, where electric field switches the sign as the 
= tt" transition even at strong coupling. This should not be taken too literally neither 
here, nor in what follows. 
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to report on the study of the continuum hmit in this framework. The surface mode picture 
[2,13] served as a reliable guide in these investigations at weak coupling. The nonpertur- 
bative nature of H-F approximation however, allows to study the model at intermediate 
and strong couplings as well. 

For the case of a single flavour, I plot in Fig. 12 the qualitative behaviour of the 
concluded full phase diagram of the model in the Hartree-Fock approximation. There is 
a planar region, embedded in the M — g'^ plane, where parity is spontaneously broken. 
From the point of view of the surface mode picture, it can be understood as the surface 
of = tt" transitions. Entering the region from the M-direction is accompanied by the 
appearance of the surface charges and the background electric field ( = — > ^ = ±7r" ) . 
Crossing the region in the Ms-direction corresponds to reversing the sign of the electric field 
("6 = TT $ = — tt"). In accordance with Aoki's scenario, tp'-f^tp acquires an expectation 
value in the broken region. Taking into account the infinite-coupling result of [17], I expect 
the parity- violating phase to extend all the way to ^r^ ^ oo. 

While the = tt <-> ^ = — tt" transitions are naturally first order, there is a strong 
evidence that the parity-violating phase transitions in M — plane are also first order in 
the Hartree-Fock approximation at nonzero g. This seems quite unnatural since one would 
have the first order ends at the boundaries of the parity-violating region. On the other 
hand, if one thinks conventionally about taking the chiral continuum limit, the first order 
transition looks quite appropriate. There are two parts to the conventional wisdom about 
taking this limit, which is usually thought about in analogy to QCD in four dimensions. 
First, since the gauge coupling constant is dimensionful (inverse length) in two dimensions, 
it is assumed that the continuum limit can only be taken at vanishing g (the dimensionless 
lattice coupling). Second, the chiral limit is assumed to be taken by following the line of 
phase transitions Mc{g'^) towards ^ = 0. The immediate consequence of these assumptions 
is that g{a) vanishes at a = and is an increasing function in the vicinity of this point. 
Here a is the lattice spacing. Denoting by A(a) the dimensionless mass gap along the 
line of phase transitions, the physical mass of the lightest particle in the theory is given 
by A (a) /a. Since chiral symmetry is broken in the continuum and the lowest mass is 
nonzero (e/y^), the above ratio should approach the constant positive value as a — > 0. 
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Consequently, similarly to g{a), A(a) should also vanish at a = and increase in the 
vicinity of this point. Puting the above two conclusions together, A.{g) must have this 
local property as well. In particular, it is zero at (7 = 0, but increases as g becomes finite. 
Hence, if the conventional picture about chiral continuum limit is correct, the transitions 
should become first order as the gauge coupling is turned on. 

In the light of the above considerations, it is not entirely obvious that the H-F ap- 
proximation is giving an incorrect answer here (which it of course well can). One possible 
solution is that what is depicted in Fig. 12 is not all that happens in the model. In par- 
ticular, there could be another sheet of first order phase transitions going off the M — g"^ 
plane and crossing this plane at Mc{g'^). If that was the case, then the second order ends 
at Mc{g^) would not appear to be necessary. However, I have not found the evidence that 
would support this scenario in the H-F approximation. In summary, the order of parity 
violating phase transition is a very interesting issue by itself. However, it can only be 
satisfactorily settled by accurate calculation beyond the H-F approximation. 

Similarly to the one-flavour case, there is a surface of = tt" transitions also in the 
model with two flavours. It takes a more complicated shape however and its qualitative 
behaviour in the H-F approximation is depicted in Fig. 13a. The "tube" of phase transitions 
encloses the (7^-axis as the explicit flavour breaking is taken to zero and touches the M — g^ 
plane at Mc{g^). Since the electric fleld generated by the surface charges switches the 
sign, the transitions across the surface of the "tube" are naturally flrst order. In the 
H-F approximation, this is so even when approaching Mc{g^) at flnite coupling. The 
transition close to Mc{g'^) becomes second order in H-F approximation only as ^ 0. 
Note also that contrary to the single-flavour case, the possible continuous nature of the 
phase transitions along M(.{g^) would not be in conflict with the above argument concerning 
the chiral continuum limit. This is because in the multi- flavour case, the nonsinglet part 
of the flavoured chiral symmetry is not anomalous and there is a massless particle in the 
continuum theory. 

The crucial difference between Fig. 12 and Fig. 13a is that with two flavours, it is only 
Mc{g^) that is shared by the M — g'^ plane and the surface of "^ = tt" transitions. If the 
identiflcation of the "M5" -physics on the lattice and the "^"-physics in the continuum is 
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correct at weak coupling, nothing special should happen upon crossing Mc{g^) with respect 
to parity, and it indeed doesn't. At strong coupling however, the above scenario might 
well break. This is nicely observed in the H-F approximation. In particular, the numerical 
evidence suggests the existence of rather strong coupling Qs (4 < (7^ < 6), so that for g > gg, 
the expectation value of ip^5Ts'ip is nonzero inside the "tube" of Fig. 13a. Consequently, 
parity-flavour is broken in M — plane at subcritical masses and strong couplings. I 
conclude the qualitative H-F phase diagram in this plane as depicted in Fig. 13b. The full 
lines in this phase diagram represent Mc{g'^) and they are also characterized by the fact 
that (V'^) exhibits a jump as they are crossed. Parity-flavour however, is only broken in 
the "BP" -region, bounded from bellow by the dashed line. 

The above results indicate that Aoki's scenario is not realized in the two-flavour 
Schwinger model at the H-F level. It is quite feasible, that this is the case for QCD in 
four dimensions as well [21,22]. Although the analogy between QED2 and QCD^ should 
certainly not be taken too seriously (especially in case of an approximation), I believe that 
the phase diagram of Fig. 13b indeed represents a possible toy picture of what might be 
happening in the latter case. In particular, that the parity-flavour broken phase shrinks 
to zero width before entering the vicinity of the continuum limit. There would still be a 
line Kc{g'^), running up from the QCD flxed point, on which a transition in (il^il^) could be 
observed. However, {ip^^Tsi/j) would remain zero. 

Appendix 

In this appendix I will briefly describe a technical detail on the numerical procedure 
used to solve the H-F equations (37). The standard way to proceed is to iteratively generate 
the sequence of sets of one particle states {S^°\ S^^\ . . . S^''\ . . so that ^('=+1) is the 
eigenset of H^'^\ Here ^^^^ is the H-F Hamiltonian with direct and exchange potentials 
determined from wavefunctions of S^'^\ With a reasonable choice of the initial set the 
sequence frequently converges well to the self-consistent set Shf- 

Note however, that there is certain ambiguity in the procedure that might cause a 
problem. It arises because of the fact that the energy of the H-F vacuum is not just a 
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sum of the one particle energies of the filled levels. Indeed, at nonzero coupling there is an 
interaction part contributing to the total energy and it may well be that the filled levels are 
not those from S^^ with lowest one-particle energies. If that is the case and the iteration 
proceeds by filling the lowest levels at each step, the procedure can never converge to a 
self-consistent set. 

This bad looking fiaw can however be quite easily rectified [19]. Instead of the original 
H-F problem (24,37), consider the one with the two body potential shifted by a constant, 
i.e. 

\ n — m \ — > \ n — m \ + C . (42) 

One naturally expects that a resulting H-F vacuum will not be physically different from 
that of the original problem. Indeed, it can be easily checked that both problems share 
their solutions. However the vacuum energy and also the one particle energies will change. 
In particular, 

E — > E + CL^ 

( 2 (43) 
^ |e"+C^, a filled; 

i e"*, a empty. 

Note that it is only the filled levels that get shifted in energy, not the empty ones. 
Therefore by choosing C to be negative and sufficiently large one can allways make the 
filled levels to be those with lowest one-particle energies. The H-F problem with such C 
can then in principle be solved by standard iteration as described above. If the solution is 
found, it is also the H-F vacuum of the original problem. 

In an actual computation, the constant C was chosen by trial and error. If the iteration 
failed for a given C, a larger value has been set. In general, larger values were needed for 
larger values of ^, as one would expect. 
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Electric Field vs Fermion Mass 
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Fig.l Vacuum expectation value of the electric field in the middle of the system 
as a function of fermion mass. 
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Electric Field vs Position 
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Fig.2a (Upper) Spatial dependence of the electric field in symmetric (M = 1.04) 
and broken (M = 0.24) phases. Fig.2b (Lower) Spatial distribution of the cor- 
responding H-F levels. 
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Fig.3 Hartree-Fock phase diagram in M — gr^ plane as seen on a finite lattice. 
Diamonds sample the critical line Mc{g^)i with parity broken in the left region. 
The square marks Mc(0) at infinite volume limit. The left vertical line represents 
Mc(oo) quoted in [17] for the model in standard Lagrangian formulation. 
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Electric Field and Axial Density vs Fermion Mass 



0.4 h 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




Electric Field « 
Axial Density + 

g2=0.1 
L = 32 



0.2 0.4 0.6 0.8 

M 



0) 
O) 

CO 

O 

"TO 



0.5 



0.4 



0.3 



0.2 



0.1 



Axial Charge Density in the Vacuum 



"T 1 1 r 



L = 32 



g^ = 1.0 



~i r 



M = 0.24 
M = 1.04 



_l I I L_ 



_l L_ 



10 15 20 25 30 
n 



Fig.4a (Upper) Electric field and axial charge density against the fermion mass. 
Transitions seem to occur simultaneously. Fig. 4b (Lower) Spatial dependence of 
axial density in symmetric (M = 1.04) and broken (M = 0.24) phases. 
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Electric Field and Axial Density in Broken Phase 
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Fig.5 Expectation values of electric field and axial charge density in broken phase. 
Relative size of the condensates varies with gauge coupling. 
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One Particle Energies vs Position 




Fig. 6 Spatial distribution of the H-F levels in broken phase at = 1.5. Filled 

levels are marked by the diamonds and empty ones by the crosses. 
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Gap vs Fermion Mass 
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Fig. 7a (Upper) Inverse correlation length from electric field and axial charge 
correlation functions. Fig.Tb (Lower) Inverse correlation length of the electric 
field correlator along the line of phase transitions. 
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Fig.Sa (Upper) Assignements of the angles dp (p) and dd (d) of Eq. (41) to a 
general point A in the M — M5 plane. Fig. 8b (Lower) Expected phase diagram 
for two flavours at weak coupling based on the surface mode picture of Ref. [2] . 
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Fig.9a (Upper) Electric field vs M5 in the two- flavour model at = 0.1 and 
M = on a lattice with 40 sites. Fig. 9b (Lower) Electric field along the M-axis 
in the same situation. 
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Fig. 10 Phase diagrams of the two- flavour model in the M — M5 plane at = 0.1 
and ^r^ = 1.2 on a lattice with 32 sites. The points represent the positions where 
electric field in the Hartree-Fock vacuum changes its sign. 
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Parity-Flavour Order Parameter vs Mass 
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Fig. 11a (Upper) Expectation value of V'Ts'T'sV' 9^ = 1-2 and M5 = 0.01 on a 
lattice with 40 sites. Fig. lib (Lower) The same expectation at g'^ = 6.0 on a 
lattice with 32 sites. 
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Fig. 12 The concluded qualitative behaviour of the full H-F phase diagram for 
the one-flavour Schwinger model on the lattice with Wilson fermions. Parity is 
spontaneously broken in the black sheet, embedded in the M — plane. 
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Fig. 13a (Upper) The concluded qualitative shape of the surface of "6* = tt" 
transitions for the model with two flavours. Fig.llb (Lower) The H-F phase 
diagram of this model in M — ^r^ plane. Parity-flavour is spontaneously broken 
in the "BP" region. 
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